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Using the BES-II detector at the BEPC Collider, we measured the lowest order cross sections 
and the R values (R — a< ^+ e -^i- l ^ Ions /°'e+e-^ f i+ f i-) f° r inclusive hadronic event production at 
the center-of-mass energies of 3.650 GeV, 3.6648 GeV and 3.773 GeV. The results lead to R u d s = 
2.224 ± 0.019 ± 0.089 which is the average of these measured at 3.650 GeV and 3.6648 GeV, and 
R = 3.793 ± 0.037 ± 0.190 at = 3.773 GeV. We determined the lowest order cross section for 
^(3770) production to be cr^ (3770) = (9.575 ±0.256±0.813) nb at 3.773 GeV, the branching fractions 
for ^(3770) decays to be BF(V>(3770) -> D°D°) = (48.9 ± 1.2 ± 3.8)%, BF(ip(3770) -> D+D') = 
(35.0±1.1±3.3)% and BF(ip(3770) -> DD) = (83.9±1.6±5.7)%, which result in the total non-DD 
branching fraction of V(3770) decay to be 5^(^(3770) -» non - DD) = (16.1 ± 1.6 ± 5.7)%. 



I. INTRODUCTION 



proceeds via quark-antiquark pair production where the 
photon couples directly to the charge of the pointlike 
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quarks. A consequence of this picture is that the to- 
tal lowest order cross section, cr^ ad , for inclusive hadron 
production in e + e~ annihilation must be proportional to 

the lowest order cross section, a® + , _ , for muon pair pro- 

fj, fi 

duction, which results in the relation 

N f 

i 

where Qi is the charge of the ith quark; the factor of 3 
accounts for three different colors of quarks; the sum runs 
over all quark flavors, Nf, involved, for which the quark 
pair production thresholds are below the e + e~ annihila- 
tion energy. The Eq. (1) indicates the ratio 

B N i 

R=^ = zJ2Q 2 i (2) 

P+M i 

to be constant as long as the cm. (center-of-mass) energy 
E cm does not overlap with resonances or the threshold of 
the production of new quark flavors. It also indicates that 
the R u ds value for continuum light hadron (containing u, 
d and s quarks) production should tend to be constant 
in the energy region above 2 GeV. 

This naive theoretical prediction for the R value has 
to be modified to take into account the finite mass of the 
quarks and the emission of the gluons by the quarks. In 
principle, the R values can be computed in the pQCD 
(perturbative QCD) with these corrections. So precise 
measurements of R values at low energy region are im- 
portant for the test of the prediction by the pQCD in 
this energy region. Moreover the R values at all energies 
are needed to calculate the effects of vacuum polarization 
on the parameters of the Standard Model. For example, 
the dominant uncertainties in the quantities a(M|), the 
QED running coupling constant evaluated at the mass of 
Z , and a M = (g — 2)/2, the anomalous magnetic moment 
of the muon, are due to the calculation of hadronic vac- 
uum polarization 1]. A large part of uncertainty in the 
calculation arises from the uncertainties in the measured 
inclusive hadronic cross sections in open charm threshold 
region, in which many resonances overlap. To get cred- 
ible measurements of R and various lowest order cross 
sections in this region, the overlapping effects have to be 
clarified clearly. 

On the other hand, the measurements of the R values 
below and above the threshold of DD production can be 
used to determine the branching fractions for -0(3770) — > 
D°D°,D+D-,DD, and for V(3770) -► non-DD with 
the measured cross sections for the D°D° and D + D~ 
together. The V(3770) resonance is believed to decay 
predominantly into DD 2]. However, there are dis- 
crepancies between the measurements of the DD cross 
section and the measurements of 0(3770) cross section 
which can be obtained from ^(3770) resonance parame- 
ters. In recent days, there are some publications to re- 
port the observation of non— DD decays of 0(3770) res- 
onance 31 41 51. In this Letter, we Drcsent more Drecise 



measurements of the R values at the cm. energies of 
3.650, 3.6648 and 3.773 GeV. With the measured R val- 
ues and the previously measured cross sections for DD 
production at 3.773 GeV we determine the branch- 
ing fractions for ^(3770) D°D°, D + D~ ,DD, and for 
V>(3770) -> non-DD 



II. BES-II DETECTOR 

The BES-II is a conventional cylindrical magnetic de- 
tector that is described in detail in Ref. pfl. A 12-layer 
Vertex Chamber (VC) surrounding the beryllium beam 
pipe provides input to the event trigger, as well as co- 
ordinate information. A forty-layer main drift chamber 
(MDC) located just outside the VC yields precise mea- 
surements of charged particle trajectories with a solid an- 
gle coverage of 85% of 4ir; it also provides ionization en- 
ergy loss (dE/dx) measurements which are used for parti- 
cle identification. Momentum resolution of 1.7%-y/I + p 2 
(p in GeV/c) and dE/dx resolution of 8.5% for Bhabha 
scattering electrons are obtained for the data taken at 
y/s = 3.773 GeV. An array of 48 scintillation counters 
surrounding the MDC measures the time of flight (TOF) 
of charged particles with a resolution of about 180 ps for 
electrons. Outside the TOF, a 12 radiation length, lead- 
gas barrel shower counter (BSC), operating in limited 
streamer mode, measures the energies of electrons and 
photons over 80% of the total solid angle with an energy 
resolution of (Je/E = 0.22/a/E (E in GeV) and spatial 
resolutions of = 7.9 mrad and oz = 2.3 cm for elec- 
trons. A solenoidal magnet outside the BSC provides a 
0.4 T magnetic field in the central tracking region of the 
detector. Three double-layer muon counters instrument 
the magnet flux return and serve to identify muons with 
momentum greater than 500 MeV/c. They cover 68% of 
the total solid angle. 



III. MEASUREMENT OF THE OBSERVED 
HADRONIC CROSS SECTIONS 

For a sample of data taken at cm. energy E cm ^, the 
observed cross section for the inclusive hadronic event 
production is obtained by 

<T° bS (D cm i) = jWggmj) (3) 

^{-CJcm.i) Chad (, -e'en!, i J e had 

where i is the ith energy point at which the data were 
collected, Nh a d(E cm ,i) is the number of the inclusive 
hadronic events observed at this energy; L(E cm ^) is the 
integrated luminosity of the data collected; ehad(D cmj i) 
is the efficiency for detection of the inclusive hadronic 
events, and e^"| is the trigger efficiency for collecting the 
hadronic events in on-line da.ta, acquisition. 
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A. Measurement of luminosity 

The integrated luminosities of the data sets are deter- 
mined by 



600 



L(E C 



N e + e -(E cm>i ) - n b 



{E c 



(^cm,t) e e +t 



(4) 



where N e + e - (E cm ^) and rib are the number of the se- 
lected Bhabha events and the number of background 
events respectively, e e + e - (E cmi ) is the efficiency for de- 
tection of the Bhabha events, e t " s _ is the trigger effi- 
ciency for collecting the Bhabha events in on-line data 
acquisition. For the data used in the analysis, the trigger 
efficiency is = (100.0±^)% (see subsection D). 

To select the candidates for Bhabha scattering e + e~ — > 
(7)e + e~, it is first required that exactly two charged 
tracks with total charge zero be well reconstructed. For 
each track, the point of the closest approach to the beam 
line must have the radius < 1.5 cm and \z\ < 15 cm 
where \z\ is measured along the beam line from the nom- 
inal beam crossing point. Furthermore, each track is re- 
quired to satisfy |cos6>| < 0.7, where is the polar an- 
gle of the charged track, to ensure that it is contained 
within the barrel region of the detector. Next, it is re- 
quired that the energy deposited for each charged track 
in BSC be greater than 1.1 GeV (i.e. ££agk >1A GeV ) 
and at least the magnitude of one charged track momen- 
tum be greater than 0.9£'b, where E^ is the beam en- 
ergy. Figure n shows the distribution of the energies de- 
posited for muons (hatched histogram) and electrons or 
positrons (points with error bars) in the BSC, where the 
data sample of the muons and the electrons or positrons 
are selected from the decays of ip(2S) — > J/i{)ir + ir~ , 
and J/ij) —* or e + e _ . From Fig. ^ one can 

see that the criterion E^^ > 1.1 GeV separates the 
e + e~ — > (7)y« + /!~ from the Bhabha scattering effectively. 
In addition, because the Monte Carlo simulation does 
not model the energy deposited well in the rib regions 
of the BSC, any charged track from the selected Bhabha 
events is required to hit one of the four regions of the 
BSC (selected z regions in BSC): (1) z B sc < -1.04 m, 
(2) -0.77 m < zbsc < -0.1 m, (3) 0.1 m < z BSC < 0.77 
m, (4) zbsc > 1-04 m. 

The two oppositely charged tracks go in opposite direc- 
tions in the R — <f> plane. Because the tracks are bent in 
the magnetic field, the positions of the two shower clus- 
ters in the R — (f> plane of the BSC are deviated from the 
back-to-back directions. We define the angle difference of 
the two clusters by 6(f) = |0i — 02 1 — 180° in degrees, where 
the 4>\ and <f>2 are the azimuthal angles of the two clusters 
in the BSC. Figure [2] shows the 6<p distribution for the 
events which satisfy the selection criteria for the Bhabha 
events. These events are from a portion of the data taken 
at 3.773 GeV. Using a double Gaussian function plus a 
second order polynomial to fit the distribution, we obtain 
the number of the candidates for e + e~ — * ("f)e + e~ . The 
acceDted candidate events are examined for backeround 
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FIG. 1: The distribution of energies deposited for muons 
(hatched histogram) and electrons or positrons (points with 
error bars) in BSC, where the sample of the muons and 
the electrons or positrons is from the decays of t/>(25) — > 
J/tpTr + TT~ , and J/if) — + or e + e~; the open histogram is 

for the electrons or positrons from the Monte Carlo events of 
ip(2S) -» J/^tt+tt-, and J/ip -> e+e~. 



contaminations by visual scan. The detailed scans for 



the accepted e H 



(j)e + e events show that about 



0.5% of the accepted events are due to background con- 
tamination. After subtracting the background, the pure 
number of e + e~ — > (^)e + e~ events is retained. 

The detection efficiency e e + e - for the Bhabha scat- 
tering e + e~ — > ("/)e + e~ is determined by analyzing the 
Monte Carlo events of e + e~ — ► (j)e + e~. These events 
are generated with the radiative Bhabha generator || 
written by R. Kleiss et at, which includes hard photon 
emission and a 3 radiative correction. 

Using the pure number selected, the visible cross sec- 
tion <r e + e - read from the generator, the detection effi- 
ciency for e + e~ — > (-f)e + e~ obtained by Monte Carlo 
simulation, and the trigger efficiency, we can determine 
the integrated luminosity of the data from Eq. (4). Ap- 
plying the procedure to the data sets taken at the three 
energy points, we get the measured integrated luminosi- 
ties of the data sets. The second column of table [I] lists 
the integrated luminosities of the data sets, where the er- 
rors are combined from statistical and systematic errors. 

The systematic uncertainty in the measured values of 
the luminosities arises mainly from the difference between 
the data and Monte Carlo simulation. Table [H] summa- 
rizes the systematic uncertainties due to the e + e~ — > 
(7)e + e~ event selection criteria. The total uncertainty in 
the measured luminositv is estimated to be about 1.8%. 
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FIG. 2: The distribution of the . 
the selected e + and e~ tracks. 
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TABLE II: The relative systematic uncertainties in the mea- 
sured luminosity due to the e + e~ — > (^y)e + e~ event selection 
criteria. 



Criterion 


% [%\ 


Radius < 1.5 cm 


0.18 


\z\ < 15 cm 


0.48 


\cos6\ < 0.7 


0.65 


Ensc > 1-1 GeV 


0.27 


P| rack (or p* rack ) > o.9E h 


1.33 


Selected z regions in BSC 


0.89 


Total uncertainty in Bhabha event 


1.83 



2.0 ns < Ttof < T p + 2.0 ns, where Ttof is 
the time-of-fiight of the charged particle, and T p is 
the expected time-of-fiight of proton with the given 
momentum; 

the charged track must not be identified as a muon; 

for the charged track, the energy deposited in BSC 
should be less than 1.0 GeV. 



TABLE I: Summary of the luminosities of the data sets, the 
numbers of the selected candidates for e + e~ — > hadrons 
and the estimated numbers of the events of the processes 
e + e~ — > l + l~ (I — T,e,n), e + e~ — » e + e~l + l~ and e + e~ — > 
e + e~ hadrons which were misidentified as the events of 
e + e~ — > hadrons. 





L 


ATZfit 

"had 


ni+i- 


n e+e-l+l- 


(GeV) 


[nb- 1 ] 






& n e+a - h 


3.650 


5537.7 ± 102.3 


54576 ± 239 


2038 


219 


3.6648 


998. 2± 19.2 


9615 ± 100 


382 


40 


3.773 


17300.0 ±319.6 


274021 ± 538 


8603 


701 



B. Selection of hadronic events 

In order to effectively remove the e + e~ — > (7)e + e _ 
and e + e~ — > ("/)fi + fi~ events from the selected hadronic 
event sample, the hadronic events are required to have 
more than 2 good charged tracks, each of which is re- 
quired to satisfy the following selection criteria: 

• the charged track must be with a good helix fit 
and the number of dE/dx hits per charged track is 
required to be greater than 14; 

• the point of the closest approach to the beam line 
must have radius < 2.0 cm; 

• | cos ^| < 0.84, where 6 is the polar angle of the 
charged track; 

• P < Eb + 0.1 x Eb x + El), where p is the 
charged track momentum and E^ is the beam en- 
erffv in GeV: 



In addition, the total energy deposited in BSC should be 
greater than 28% of the beam energy. Furthermore, the 
selected tracks must not all point into the same hemi- 
sphere in the z direction. No criterion for the number of 
the observed photons is applied to the selected hadronic 
events. 

Some beam-gas associated background events can also 
satisfy above selection criteria. However, the beam-gas 
associated background events are produced at random z 
positions, while the hadronic events are produced around 
2 = 0. This characteristic can be used to distinguish 
the hadronic events from the beam-gas associated back- 
ground events. To this end, the averaged z of the charged 
tracks in each event is calculated. Figure shows the 
distribution of the averaged z. These events are from 
a portion of the data taken at 3.773 GeV. In Fig. [31 
the points with error bars show the events from the 
Monte Carlo sample which is generated with the genera- 
tor described in Section III.C and simulated with the 
GEANT3-based Monte Carlo package , the histogram 
shows the events from the data, and the shadowed his- 
togram shows the events from the separated beam data. 
Using a Gaussian function plus a second order polynomial 
to fit the averaged z distribution of the events, we obtain 
the number, iVj^, of the candidates for hadronic events. 
The third column of table [I] lists obtained from the 
data sets taken at each of the energy points, where the er- 
rors are combined from statistical and systematic errors. 
This number of candidates contains some contaminations 
from some background events such as e + e~ — > t + t~, 
e + e~ — > (7)e + e~, e + e~ — ► (7)/i + ^ - and two-photon ex- 
change processes. The number of the background events 
can be estimated by means of Monte Carlo simulation 
(see Section TTT.FA 
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TABLE III: The relative systematic uncertainties in measur- 
ing the produced hadronic events due to the event selection. 



1 Lei 1UI1 


a i //v prd r%i 

£i Jv prd/JV had [/Oj 
had 


Good helix fit 


0.10 


N dE/dx > 14 


0.01 


K H < 0.02 m 


1.30 


\cos6\ < 0.84 


0.36 


p < E b + 0.1 x E b x 7(1 + S 6 2 ) 


0.07 


2.0 rzs < T TO f < Tp + 2.0 ns 


0.43 


/xld = .false. 


0.04 


SbIc < 1.0 GeV 


0.56 


Ebsc > 0.28£ b 


0.49 


Same hemisphere cut 


1.01 


Fitting to averaged Z distribution 


0.65 


-Wchtrk > 3 


1.5 


Total uncertainty 


2.50 




-25 25 50 75 

Averaged Z [cm] 



FIG. 3: The distribution of the averaged z of the charged 
tracks which satisfy the hadronic event selection criteria, 
where points with error bars show the events from the Monte 
Carlo sample, the histogram shows the events from the data, 
and the shadowed histogram shows the events from the sepa- 
rated beam data; the curves give the best fit to the z distri- 
bution. 



C. Monte Carlo method and the efficiency thad 

Due to ISR (Initial State Radiation), the effective cm. 
energy for the e + e~ annihilation is E c ff — J s(l — x), 
where xy r s/2 is the total energy of the emitted photons 
and yfs is the nominal cm. energy. For a certain en- 
ergy point in experiment, the experimentally observed 
hadronic events are not only produced at the i/s, but 
produced in the full energy range from the y/s to ~ 0.28 
GeV (for production of two pions). To determine the ef- 
ficiency for detection of hadronic events produced in the 
full energy range, we developed a special Monte Carlo 
generator in which the initial state radiative correc- 
tion to a 2 order is taken into account. 

Figure 0] shows the differential cross section d(r/dE e g 
for the inclusive hadronic event production when set- 
ting the nominal cm. energy to be at 3.80 GeV. At 
an effective cm. energy, the final hadronic states are 
produced by calling the sub-generators such as LUND 
model and the resonance generators including 

^(3770), ip(2S) } J/tp E3, 0(1020), p(770), and u(782) 
etc. according to the corresponding lowest order cross 
sections of these processes, respectively. The CMD-2 
7T7T production data [l4| with Gounaris-Sakurai param- 
eterization are used to simulate the spectrum of 
p(770) and the p-u> mixing in the energy range below 
1.2 GeV. The resonances are set to decay into all possi- 
ble final states according to the known decay modes and 
branching fractions. These generated events are simu- 
lated with the GEANT3-based Monte Carlo simulation 
package. The reconstructed Monte Carlo events are then 
fed into the analysis program to determine the efficien- 
cies, Chad, for measurements of the observed cross sec- 
tions for inclusive hadronic event production at each of 
the three energy points. 

For simulations of the inclusive hadron production, pa- 
rameters in the LUND generator are tuned using an in- 
clusive hadronic event sample of 5.5 x 10 5 events from 
the data taken at 3.65 GeV with the BES-II detector. 
The parameters are adjusted to reproduce good agree- 
able distributions of some main kinematic variables be- 
tween data and Monte Carlo sample. The uncertainty in 
£had due to the adjusted parameters is estimated to be 
~ 0.6%, while the uncertainties due to the errors of the 
•0(3770) and tp(2S) resonance parameters are estimated 
to be ~ 1.5% and ~ 1.2%, respectively. Combining these 
uncertainties in quadrature yields the systematic uncer- 
tainty in the efficiency e^ad to be about 2%. 

The second column of table HVl lists the efficiencies for 
detection of the inclusive hadronic events at three energy 
points in the case of setting the continuum R value to be 
at 2.26 (see section IV). 



The systematic uncertainty in measuring the produced 
hadronic events due to the hadronic event selection crite- 
ria is estimated to be about 2.5%. Table ITTT1 summarizes 
the relative systematic uncertainties in selecting the pro- 
duced hadronic events. 



D. Trigger efficiency 

The requirements of the trigger for recording the data 
on-line axe almost the same as those used in collecting 1 the 
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TABLE IV: Summary of the efficiencies for detection of the 
inclusive hadronic events and the observed cross sections for 
e + e~ — ► hadrons at three energy points, where the error in 
Qiad is statistical only; while the errors in o^ad are statistical 
and point-to-point systematic, respectively. 





£had 


obs 
°had 


(GeV) 




(nb) 


3.650 


0.4977 ± 0.0022 


18.983 ±0.087 ±0.182 


3.6648 


0.5033 ± 0.0023 


18.299 ±0.199 ±0.176 


3.773 


0.5528 ± 0.0022 


27.680 ± 0.056 ± 0.266 



Nh&d{E cm ,i) with ^ad _ N b, where Nb is the number of 
the background events, such as e + e~ — > t + t~ , e + e~ — > 
(7)e + e~, e + e~ — > (7)/i + /i~ and two-photon exchange 
processes. The number of the background events can be 
estimated by using the theoretical cross sections of these 
processes, the rates of misidentifying these processes as 
hadronic events and the total integrated luminosities of 
the data sets, which is given by 

N b = L X (yii+i-(Ti+i- +Ve+e-l+l-°'e+e-l+l- 
~l - ^7e+ e~ 7T+7T - *^e+ e~ 7r+7r — 

~l _ ^?e+ e~ hadrons® e~ hadrons) (^) 




0.5 1 1.5 2 2.5 3 3.5 4 

E eS [GeV] 

FIG. 4: The differential cross section for the inclusive 
hadronic event production when setting the nominal cm. en- 
ergy to be at 3.80 GeV; the histogram shows the resonances 
and continuum hadronic event production; the dashed line 
shows the cross sections for fi + fj,~ pair production. 

data for the work and the work However, for the 
-0(3770) data acquisition, we slightly modified the trig- 
ger requirements for the charged tracks, which results in 
a little bit improvement in recording the two charged trak 
events. The trigger efficiencies are obtained by compar- 
ing the responses to different trigger requirements in the 
data taken at 3.097 GeV during the time period taking 
the data at y/s = 3.773 GeV. The trigger efficiencies are 
measured to be 100.0% for both the e + e~ — > ( / f)e + e~ and 
e + e~ — > hadrons events, with an uncertainty of '\%. 

E. Observed hadronic event cross sections 

The observed cross section for inclusive hadronic event 
nroduction can be obtained from En. f3l substituting 



where &i+i-, a e + e -i+i-, cr e +e-7r+7r- and a e + e - hadrons are 
the cross sections for e + e~ — > , e + e~ — > e + e~Z + Z~, 
e + e~ — * e + e~7r + 7r~, and e + e~ — + e + e~ hadrons pro- 
cesses, respectively; while r?;+;-, n e+e - l+l -, r] e + e -„+„- , 
and r/ e + e ~ hadrons are the corresponding misidentification 
rates. 

In the calculation of the t + t~ cross section, we con- 
sider the contributions from ip(2S) decay, the QED pro- 
duction and their interference; we also consider the ef- 
fects of the initial and final state radiative corrections 
and Coulomb interaction on the cross section 0]. For 
e + e~ — > (j)e + e~ e + e~ — » ( r y)u + H~ and two- 
photon processes e + e~ —> e + e~l + l~ 20], the cross sec- 
tions are read from respective generator outputs. As 
for the estimate of the total cross section for e + e~ — > 
e + e~hadrons, we employ the equivalent photon approx- 
imation formalisms to deal with the 7-7 collision sub- 
process [2lJ H2- In the sub-process, the energy depen- 
dence of the total hadronic cross section can be described 
well by the formula of Donnachie-Landshoff parameteri- 
zation j2^| above the three pion threshold. For the con- 
tribution in the low energy region below the three pion 
threshold, it is good enough to use the simple point-like 
7r + 7r~ production cross section in the calculation [21|. 

The rates of misidentifying the above processes as the 
hadronic events are obtained from Monte Carlo simu- 
lation with the generators mentioned above. For the 
two-photon process e + e~ — > e + e~ hadrons, we use the 
Monte Cairo generator described in [2(j to simulate the 
process and determine the rates of misidentifying the 
processes as the hadronic events. The fourth and fifth 
columns of Table [fl give the estimated numbers {ni+i-, 
n e + e ~i+i- and n e + e - h ) of the background events from the 
e + e~ — ► and the two-photon exchange processes, 

which are misidentified as the inclusive hadronic events. 

Inserting N b , L, e ha d and e£j§ in Eq. (3), 

we obtain the observed cross sections for the inclusive 
hadronic event production at each of the three energy 
points, which are summarized in table ITVl where the first 
error is statistical and second point-to-point systematic 
error arising from the uncertainty in e e + e - (0.6%), un- 
certainty in N e + e - (~ (0.2 ~ 0.5%)) and uncertainty in 
£had (0.5%). The common systematic uncertainty is not 
included vet. 
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IV. LOWEST ORDER CROSS SECTION 



Radiative corrections 



with £ = -^r 2 -, where mi is the lepton mass. 

In the structure function approach by Kuraev and 
Fadin H, 



To get the lowest order cross section for the inclusive 
hadronic event production in e + e~ annihilation, the ob- 
served cross section has to be corrected for the radiative 
effects including the initial state radiative corrections and 
the vacuum polarization corrections. The correction fac- 
tor, (1 + 6(a)), is given by 



F(x,s) = [3x l3 ~ 1 5 v+s + S H , (14) 
where (3 is the electron equivalent radiator thickness, 







2a f a 



1 



(15) 



(1 + 6(8)) 



cxp 



(s) 



cr B (s) ' 



(6) 



where <r cxp (s) is the expected cross section and <J B (s) 
lowest order cross section for the inclusive hadronic event 
production. 

The expected cross section for hadronic event produc- 
tion can be written as 



rV+S 



i n a fir 2 



H 



24 



im4 



37 



(16) 



(17) 



cxp 



(s) 



dx ■ 



<x B (s(l-aQ) 
\l-U(a(l-x))\- 



;F(X,S), 



(7) 

where c B (s(l — x)) is the total lowest order cross section 
in the energy range from 0.28 GeV to yfs (or from 3.729 
GeV to y/s in the case of considering the DD produc- 
tion), F(x, s) is a sampling function and \i^yi(s}i-x))\ 2 * s 



the correction factor for the effects of vacuum polariza- 
tion including both the leptonic and hadronic terms in 
QED H3, with 



n(s ) = n had ( s ) + n ; ( s ), 



(8) 



the effects of hadronic vacuum polarization can be calcu- 
lated via the dispersion integral plj 



IIhad(s ) 



while 



47T 2 C 



]> V r ^4ds" - ina B (a) 

4m: 



„2 s — s 



Hi(s 



ls i :s(a 



2n 



/(£), 



(9) 
(10) 

(11) 



*? = -/>(!-!) 



(18) 



°2 = 



4(2 - x) In- - 



1 + 3(1 -a;) 2 



ln(l - x) - 6 



(19) 

In above expressions, m e is the mass of electron and a is 
the fine structure constant. 

For the resonances, such as ip(2S), J/tp, 4> an d w j we 
use the Breit-Wigner formula 



o*(s) 



i2Av h 



(a' - M 2 ) 2 + T 2 M 2 



(20) 



to calculate the lowest order cross section, where M and 
r are the mass and the total width of the resonance, and 
r ee = r^Jl — II(s )| -2 and Th are the partial widths 
to the e + e~ channel and to the inclusive hadronic final 
state, respectively. For the ■0(3770) resonance, we use 



i27rro e r tot ( s ) 



(s - M 2 ) 2 + T 2 ot (s')M 2 



(21) 



to calculate the lowest order cross section, where r to t (s ) 
is chosen to be energy dependent and normalized to the 
total width Ttot at the peak of the resonance. The 
rtot(s ) is defined as 



r to t(s ) = r D o£>o(s )+T D+D -(s )+T Ti 



i-DD 



(a ), (22) 



JV ° J 9 3 6 



/(0 = - 5 -| 



5 e , vr^i(2+o 



i + yr^e 
i-vr^e 



. (e < i), 

(12) 



tan 



«>1), 



rial 



where T D of,o(a ), T D + D -(s ) and T noa _ D g(s ) are the 
partial widths for ^(3770_) -> D°D°, 0(3770) -> D+D~ 
and 0(3770) — > non— DD, respectively, which are taken 
in the form 



r.DO£)o(s ) = r 8(e c 



om ^ (PD ) 3 ^ + (rp DO o) 2 
2M D o) n vq - ■ — — £>oo, 



(p° )3 l + (rp D o) 2 



(23) 
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r D +D-(s ) = r e{E c 



2Mr 



( PD+ fl + {rp° D+ f 
V D+ ) 3 1 + 0"P D+ ) S 



and 



non — DD 



0) = r [i - Boo - , 



(24) 



(25) 



where p° D and p^, are the momenta of the D mesons pro- 
duced at the peak of 0>(377O) and at the cm. energy 
a/s 7 , respectively; To is the total width of the 0(3770) at 
the peak, and r is the interaction radius of the cc, which 
is set to be 1.0 fm; Boo and are the branching frac- 
tions for ^(3770) -> D°D° and 0(3770) -> D+D~, re- 
spectively; 9(E cm — 2Mn°) and 0(E cm — 2M D +) are the 
step functions to account for the thresholds of the DD 
production. 

In the calculation of the lowest order cross section, the 
0(3770) resonance parameters M = 3772.3 ± 1.0 MeV, 
r = 25.5 ± 3.1 MeV and T ee = 0.224 ± 0.031 keV 
measured by BES Collaboration [3(j are used. Insert- 
ing the resonance parameters of J/0, ij) (2S ) quoted from 
PDG 26] and the R = 2.26 ± 0.14 [||Jl6| for the light 
hadron production in the energy range from 2.0 to 3.0 
GeV measured by BES Collaboration in Eqs. (6)-(25), 
we obtain the radiative correction factors at the three en- 
ergy points, which are summarized in tabled In deter- 
mination of the value of (1 + S(s)), the input of R value 
in calculating the cross section for continuum hadronic 
event production affects the value of (1 + S(s)). Vary- 
ing the input R by ±10% causes a variation of ±1.3% 
in (1 + 5(s)), which results in a variation of the product 
e?iad(l + S(s)) by only ±0.4%. The most uncontrolled 
cross sections in the calculation of (1 + S(s)) come from 
the hadronic cross sections in the energy range from 1.2 
to 2.0 GeV. However, the whole contribution of the cross 
section from this energy range is less than 5% of the total 
observed cross section in our case. Since the efficiency is 
quite low for detection of the hadronic events from this 
energy range (ehad < 10%), the amount of the product 
£had(l + S(s)) would also be rather stable with error less 
than 0.4%. Taking into account the uncertainty in the 
measured hadronic event production and the errors of 
the resonance parameters together, the total uncertainty 
in R measurement due to the calculation of (1 + S(s)) is 
then estimated to be less than 1.5% in this work. 



B. Lowest order cross sections and R values 

The lowest order cross section for inclusive hadronic 
event production is obtained by 



(26) 



where 



°h( s ) = a e+e-^h a drons( s ) + H °Re 8 ,i( s )> ( 27 ) 



in which ^+ e -^i iac i rons ( s ) is the cross section for inclu- 
sive hadronic event production through one photon an- 
nihilation, cr|| cs i(s) is the cross section for the ith reso- 
nance, such as J/0, ip(2S), ■0(3770) etc. which decays 
into hadronic final states. 

To obtain the lowest order cross section erP(s) 



Ruds 



<7~+ _ (s) 



^hadrons 



(s) for the hadronic 



event production through one photon annihilation at the 
energies of 3.650 and 3.6648 GeV, and the lowest or- 
der cross section crg(s) = (R u ds + #v>(3770)) ' rf+u- = 



7 b 

e+e — —thadrons 



(s) 



B 

V>(3770) 



(s) for both one photon an- 



nihilation and 0(3770) production at 3.773 GeV, the 
amount of the cross section due to the resonance pro- 
duction at the energies of 3.650 and 3.6648 GeV, and the 
amount of the cross section due to the resonance pro- 
duction but '0(3770) at 3.773 GeV have to be subtracted 
out. The third column of table |Vl summarizes the lowest 
order cross sections <J®(s). Dividing the ct°(s) by the low- 
est order cross section for /i + pT production at the same 
cm. energy, we obtain the R values, which are summa- 
rized in the fourth column of the table. The first error in 
the measured lowest order cross section and the R value 
listed in table [V] is statistical, the second is the point- 
to-point systematic and the third is common systematic 
error. 

The common systematic error arises from the uncer- 
tainty in luminosity (~ 1.8%), in selection of hadronic 
event (~ 2.5%), in Monte Carlo Modeling (~ 2.0%), in 
radiative correction (~ 1.5%) for the measured cross sec- 
tions and R values at the three energy points, and the 
uncertainty in -0(3770) resonance parameters (~ 2.7%) 
for those at 3.773 GeV only. Adding these uncertainties 
in quadrature yields the total systematic uncertainties to 
be ~ 4.0% and ~ 4.9% for the measured hadronic cross 
sections and R values for the data taken below the DD 
threshold and at 3.773 GeV, respectively. 

Averaging the R values measured at the first two en- 
ergy points (3.650 and 3.6648 GeV) by weighting the 
combined statistical and point-to-point systematic errors, 
we obtain 

R uds = 2.224 ± 0.019 ± 0.089, 

where the first error is combined from statistical and 
point-to-point systematic errors, and the second is com- 
mon systematic. This R u ds excludes the contribution 
from resonances and reflects the lowest order cross section 
for the inclusive light hadronic event production through 
one photon annihilation of e + e~. So it can be directly 
compared with those calculated based on the pQCD. The 
value is consistent with R u d s — 2.26 ± 0.14 obtained by 
fitting those |12| measured in the energy region between 
2.0 and 3.0 GeV @ and with R uds = 2.21 ±0.13 obtained 
from fitting to the inclusive hadronic cross sections for 
both the ^/>(25 l ) and 003770) resonances in the energy 
region from 3.666 to 3.897 GeV 30 . 

Figure [3] displays the values of R from this measure- 
ment and Drevious measurements bv BES Collabora- 
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TABLE V: Summary of the radiative correction factors, the 
lowest order cross sections and the R values measured at three 
energy points. 

Energy (1 + 5(a)) v°(s) R 

(O-V) (nb) 

3.6500 1.287 14.624 ± 0.067 ± 0. 142 ± 0.590 2 .243 ± 0.010 ± 0.022 ± .090 
3.6648 1.261 14. 151 ± 0. 158 ± 0. 139 ± 0.580 2 . 188 ± 0.024 ± 0.022 ± 0.088 
3.7730 1.195 23.142 ± 0.047 ± 0.222 ± 1.158 3 . 793 ± 0.008 ± 0.036 ± 0. 190 



V. BRANCHING FRACTIONS FOR THE 
DECAYS 0(3770) -> D°D° ,D+ D~ ,_DD AND FOR 
0(3770) -> non - DD 

Assuming that there are no other new structure and 
effects except the 0(3770) resonance and the continuum 
hadron production in the energy region from 3.70 GeV 
to 3.86 GeV, the branching fraction for 0(3770) — > DD 
can be determined by 



CD 



5 
4.5 

4 
3.5 

3 
2.5 

2 
1.5 



i i i 

A Gamma2 

T MARK-I 

O PLUTO 

• BES-II (PRL 84(2000)594) 

O BES-II (PRL 88(2002)1018021-1) 

★ BES-II (this work) 



T T 




3.25 3.5 

E cm [GeV] 



3.75 



FIG. 5: The values of R measured by BES Collaboration, 
MARK-I Collaboration, 77 2 Collaboration and PLUTO Col- 
laboration in the energy region between 2.85 and 3.90 GeV, 
where the error bars show the combined statistical and sys- 
tematic errors in quadrature. 



tion MARK-I Collaboration E3, 772 Collabo- 

ration |28| and PLUTO Collaboration J2<| in the energy 
region between 2.85 and 3.90 GeV. The error bars shown 
in the figure are obtained by combining statistical and 
systematic errors in quadrature. 

Using the measured R value at 3.773 GeV listed in 
table and the R u ds value for light hadron production 

V>(3770) 



measured below the DD threshold we obtain the R 
due to ■0(3770) decays to be 



R 



i/>(3770) 



= 1.569 ±0.042 ±0.133, 



(28) 



where the first error is combined from statistical and 
point-to-point systematic error and the second common 
systematic. In estimation of the systematic uncertainty, 
we assumed that the same amount of the systematic un- 
certainties in the measured values of the R and the R u ds 
is canceled in subtracting the R u d s from the R. The 
corresponding lowest order cross section for ■0(3770) pro- 
duction is 



£F(0(377O) -> DD) 



obs 
DD 



(1 + 5) 



DD 



a 



^(3770) 



(30) 



where <r° b ? and ^ (3770) 



are the observed and lowest 
order production cross sections for DD and inclusive 
hadronic events, respectively; (1 + 6) DD - is the radia- 
tive correction factor for DD production. Inserting the 
003770) resonance parameters (M = 3772.3 ± 1.0 MeV; 
T tot = 25.5 ± 3.1 MeV and T ee = 0.224 ± 0.031 keV) 
measured by BES Collaboration [3(| in Eqs. (6) and (7) 
with combining the Eqs. (7)-(25) together, we obtain the 
radiative correction factor 



(1 + *) 



DD 



0.764 ±0.014, 



(31) 



where the error is the uncertainty arising from the er- 
rors of the 0(3770) resonance parameters, the uncer- 
tainty in vacuum polarization correction and the uncer- 
tainty arising from varying the branching fraction for 
0(3770) -> DD from 84% to 100%. 

BES Collaboration measured the observed cross sec- 
tions for D°D° and D + D~ production at cm. energy 
Vs = 3.773 GeV to be a D o B o = (3.58 ± 0.09 ± 0.31) nb 
fmd<j D+D - = (2.56 ± 0.08 ± 0.26) nb @. These observed 
cross sections were obtained by analyzing the same data 
set from which the R value at y/s = 3.773 GeV is mea- 
sured. 

Inserting the 0^(3770) > the observed cross sections for 
D°D°, D+D-, DD production at 3.773 GeV and the 
radiative correction factor, (1 + S) d q, in Eq. (30), we 
obtain the branching fractions for the decays 0(3770) — > 
D°D°, D+D-, DD to be 

BF(0(377O) -> D°D°) = (48.9 ± 1.2 ± 3.8)%, (32) 



BF(0(377O) -» D+D-) = (35.0 ± 1.1 ± 3.3)%, (33) 



and 



£F(0(377O) -► DD) = (83.9 ± 1.6 ± 5.7)%, (34) 

which results in the non-DD branching fraction of 
0(3770) to be 



B 

V>(3770) 



(9.575 ± 0.256 ± 0.813) nb. (29) BF(ib(3770) -> non - DD) = (16.1 ± 1.6 ± 5.7)%. (35) 
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TABLE VI: Sources of the uncanceled systematic uncertain- 
ties in the measured cross sections for DD and -0(3770) pro- 
duction, where upper N and upper c mean the uncertainties 
for the neutral mode D°D° and the charged mode D + D~ , 
respectively; for the uncertainties in the measured cross sec- 
tion for DD production, please refer to the reference 0; we 
here re-estimate the uncertainty in the combined tracking and 
kinematic fit in selection of D events to be ~ 4.0%. 



Source 




Source 


Ao >(3770) 


[%] 


CT iK3770) 
[%] 


Particle ID 


~ 1.5 


Monte Carlo Modeling 


- 2.0 


Tracking & K.F. 


~ 4.0 


Hadron selection 


~ 2.5 


F. P. 


~ 3.0 






MC Statistics 


~ 0.6 






Br for D° 


~ 3.2 






(Br for D+) 


(-6.5) 






Total 


-6.2 N 
(~ 8.4) c 


Total 


~ 3.2 



where the first error is statistical and the second sys- 
tematic arising from uncanceled systematic uncertain- 
ties. The uncanceled relative systematic uncertainties 
are ~ 6.2%, - 8.4% and ~ 3.2% for the a^ Do , a°^ D - 
and the 0^(3770) i respectively. The systematic error also 
includes the common uncertainty of — 2.7% arising from 
the statistical uncertainty in the measured lowest order 
cross section for 0(3770) production and the uncertainty 
(~ 1.8%) in radiative correction factor (1 + S) d q. Ta- 
ble ^T] summarizes sources of the uncanceled systematic 
uncertainties for the measured (T^k and 0W377O) ■ The 
uncertainties in luminosity (~ 1.8%), in -0(3770) reso- 
nance parameters (~ 2.7%) and in radiative correction 
(~ 1.5%) are canceled out in the estimation of the sys- 
tematic uncertainty in the measured branching fractions. 



VI. DISCUSSION ABOUT INTERFERENCE 
EFFECTS 

The measured R values discussed in above sections are 
obtained based on the same treatment on the measure- 
ments of inclusive hadronic cross sections in which no 
interference between the inclusive hadronic final states 
of the resonance decays and the inclusive hadronic final 
states from non-resonance annihilation of e + e~ is taken 
into account Q^H^HH- However, since the cm. ener- 
gies of 3.650 GeV, 3.6648 GeV and 3.773 GeV are close 
to the ip(2S) resonance, there may be interference effects 
between the final hadronic states from the ip^iS) electro- 
magnetic decays and the continuum hadron production 
in annihilation of e + e~. These interference effects dis- 



TABLE VII: Summary of the radiative correction factors, 
the lowest order cross sections and the R values measured at 
three energy points, where the interference effects between the 
0(2$) electromagnetic decays and hadron production through 



non-resonant annihilation of e 



are taken into account. 



Energy 


(1 + «(«)) 


-S crr <-) 




R 


(GoV) 




(nb) 






3.6500 


1.300 


14.855 ± 0.067 ± 0.140 ± 0.584 


2 


279 ± 0.010 ± 0.022 ± 0.089 


3.6648 


1.285 


14.543 ± 0.155 ± 0.137 ± 0.570 


2 


249 ± 0.024 ± 0.021 ± 0.086 


3.7730 


1.198 


22.951 ± 0.047 ± 0.222 ± 1.155 


3 


762 ± 0.008 ± 0.036 ± 0.189 



tort the line shape of the continuum hadron production 
cross section around the ip(2S) peak. With the definition 
of the R given in Eq. (2), we can estimate the destruc- 
tive/constructive amount of the cross section due to the 
interference effects, which is given by 



interf 
'had 



R a 



interf 



(36) 



The destructive/constructive amounts of the cross sec- 
tions given in R are estimated to be —0.0581, —0.1026 
and +0.0219 at 3.650 GeV, 3.6648 GeV and 3.773 GeV, 
respectively. After correcting the cross section cr°(s) for 
the destructive/constructive amounts due to the interfer- 
ence effects, we obtained the lowest order cross section 
crr (s). The third column of table IVIll summarizes the 
a h crr ( s )- I n the case of considering the interference ef- 
fects, the correction factor (l + <5(s)) is also changed. The 
second column of tabic IVIll lists the correction factor at 
the cm. energies. Dividing the crfj crr (s) by the lowest or- 
der cross section for production at the same cm. 
energy, we obtain the R values, which are summarized 
in the fourth column of the table. The errors are sta- 
tistical, the point-to-point systematic and the common 
systematic as discussed before. 

The weighted average of the R values measured at the 
first two energy points is 

R uds = 2.268 ±0.019 ±0.091, 

where the first error is combined from statistical and 
point-to-point systematic errors, and the second is com- 
mon systematic. 

Following the same procedure as discussed in Sections 
IV and Section V, we obtained the lowest order cross 
section for 0(3770) production to be 

cr^ ( 377 0) = (9.113 ± 0.255 ±0.805) nb, (37) 

and the branching fractions for -0(3770) decays to be 

BF(^(3770) -> D Q D°) = (51.4 ± 1.3 ± 4.0)%, (38) 



BF(0(377O) -► D+D-) = (36.8 ± 1.1 ± 3.5)%, (39) 



BF(^(3770) -► DD) = (88.2 ± 1.7 ± 6.0)%. 



(40) 
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and 

BF(^(3770) -> non -£>£>) = (11.8 ± 1.7 ± 6.0)%, (41) 

where the errors are statistical and the systematic arising 
from some uncanceled systematic uncertainties. 

VII. SUMMARY 

From the analysis of the data taken at 3.650, 3.6648 
and 3.773 GeV in e + e _ annihilation, we measured the R 
values at the energy points, which are R = 3.793±0.037± 
0.190 at 3.773 GeV and R uds = 2.224 ± 0.019 ± 0.089 at 
y/s below the DD threshold. Based on the measured 
R values, we determined the lowest order cross section 
for -0(3770) production at 3.773 GeV to be 0^( 377O ) = 
(9.575 ± 0.256 ± 0.813) nb, the branching fractions for 
0(3770) decays to be BF(^(3770) -> D°D°) = (48.9 ± 
1.2±3.8)%, BF(V>(3770) -> D+D~) = (35.0±1.1±3.3)% 



and BF(0(377O) -> DD) = (83.9 ± 1.6 ± 5.7)%, which 
result in the non-DD branching fraction of ■0(3770) to be 
BF(V>(3770) -mon -£>£>) = (16.1 ± 1.6 ± 5.7)%. These 
branching fractions for 0(3770) decays are measured for 
the first time. 
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